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We used next-generation sequencing to characterise the complete mitochondrial genome of the damselfly 
Mnais costalis (Odonata, Calopterygidae). Illumina paired end reads were mapped against COI and 16S 
sequences from M. costalis and then extended using an iterative de novo map procedure. The final assem- 
bly was a contiguous sequence of 15,487 bp, which contained all standard mitochondrial coding regions 
and the putative A + T rich region. The gene configuration of the M. costalis mitogenome is similar to that 
of other odonates, comprising 13 protein-coding genes, large and small rRNA genes, and 22 tRNA genes. 
We found three intergenic spacers that are also present in all available whole odonate mitogenomes. Base 
composition of the M. costalis mitogenome is 40% (A), 20% (C), 14% (G) and 26% (T), with a high 
A +T content (66%). The characterisation of the complete mitochondrial genome of M. costalis adds to 
the growing list of mitogenomes currently available for odonates, and will help to improve primer design 
for future population genetic studies. A phylogenetic analysis including the currently available mitochon- 
drial genome sequences of odonates suggests that Epiophlebia superstes is more closely related to the 
Zygoptera than to the Anisoptera. 
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Introduction 


Despite the problems that may arise from the use of mitochondrial markers in instances of intro- 
gression or hybridisation (e.g. Hayashi, Dobata, & Futahasi, 2005), mitochondrial sequence 
data continue to provide the majority of the genetic markers used in phylogenetic studies on 
odonates (Ballare & Ware, 2011). This widespread use of mitochondrial markers for such stud- 
ies is mainly due to the limited number of alternative nuclear genes that could be amplified 
with nearly universal primers across different species (but see Ferreira et al., 2014). Having 
full mitochondrial genomes available provides the opportunity to improve primer design, and 
also to increase the number of target genes used in phylogenetic inference, therefore providing 
more robust phylogenetic reconstructions (e.g. Lin, Chen, & Huang, 2010). Also, mitochondrial 
genomes may allow diverse comparative and evolutionary genomics questions to be answered, 
not only in odonates, but in insects in general; such as the evolution of genome size, the study of 
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Table 1. Mitochondrial genomes currently available for Odonata. 





Mitogenome GenBank 
Suborder Family Species size (bp) Accession Reference 
Anisoptera Libellulidae Orthetrum triangulare 14,033* AB126005 Yamauchi, Miya, M. U., 
melania & Nishida, 2004 
Hydrobasileus 15,088 KM244659 Tang et al., 2014 
croceus 
Brachythemis 15,056 KM658172 Yu, Cheng, Ma, Yu, & 
contaminata Zhang, 2016 
Corduliidae Cordulia aenea 14,448* JX963627 Simon & Hadrys, 2013 
Gomphidae Davidius lunatus 15,913 EU591677 Lee et al., 2009 
Ictinogomphus sp. 15,393 KM244673 Tang et al., 2014 
Anisozygoptera Epiophlebidae Epiophlebia superstes 15,435 JX050223 Wang et al., 2015 
Zygoptera Platycnemidae Platycnemis foliacea 15,382 NC027180 Unpublished 
Coenagrionidae Ischnura pumilio 15,250 KC878732 Lorenzo-Carballa et al., 
2014 
Euphaeidae Euphaea formosa 15,700 HM126547 Lin et al. 2010 
E. yayeyamana 15,709 KF718293 Unpublished 
E. ornata 15,863 KF718295 Unpublished 
E. decorata 15,861 KF718294 Unpublished 
Pseudolestidae Pseudolestes mirabilis 15,122 FJ606784 Unpublished 
Calopterygidae Atrocalopteryx atrata 15,424 KP233805 Unpublished 
Vestalis melania 16,685 JX050224 Chen et al., 2015 


*indicates that the mitogenome sequence is incomplete. 


gene rearrangements (Cameron, 2014), or to analyse the evolutionary processes that underlie the 
evolution of synonymous and non-synonymous codon usage (e.g. Kumar et al., 2012). 

In the last couple of years, there has been a notable increase in the number of mitochondrial 
genome data available for odonates; and currently the complete (or nearly complete) mitochon- 
drial genomes of 16 odonate species belonging to nine families are available (Table 1). However, 
and despite it covering nearly a third of the families currently described within the order, this 
constitutes yet a small fraction of the extant odonate diversity (~ 6000 species in 30 families; 
Dijkstra et al., 2013). 

Here, we report the complete mitochondrial genome of the damselfly Mnais costalis 
(Zygoptera, Calopterygidae) from Japan, obtained by next generation sequencing. Members of 
the Calopterygidae, commonly known as broad-winged damselflies, demoiselles or jewelwings, 
have metallic-coloured bodies and some species also have conspicuously pigmented wings. 
They have been used as models for research into reproductive behaviour, interspecific inter- 
actions, character displacement and sperm competition (Cérdoba-Aguilar, 2008), with Mnais 
damselflies being particularly important as they are one of the few species that have a male-linked 
colour polymorphism that is associated with a behavioural phenotype (Plaistow & Tsubaki, 
2000; Tsubaki, Hooper, & Siva-Jothy, 1997; Tsubaki & Okuyama, 2016). Finally, we reconstruct 
the phylogenetic relationships among odonates, based on the mitochondrial genomes currently 
available. 


Material and methods 


Sample collection, DNA extraction and sequencing 


Five adult Mnais costalis males (two orange and three clear winged) were collected at Togichi 
Prefecture (36°42’ N, 140°13’ E) and used to extract genomic DNA from the thoracic muscle. For 
DNA extraction, we used the Qiagen Genomic-Tip 500/G (Venlo, The Netherlands), following 
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the manufacturer’s instructions. Samples were pooled for sequencing (pool 1 — clear and pool 2 
— orange) on a HiSeq 2000 (Illumina, California, USA), generating 2 x 100 bp paired end reads 
(approximate fragment size 550-750 bp), with chemistry v.3, at the Centre for Genomic Research 
(www.liverpool.ac.uk/genomic-research/). Information on the number of reads obtained and data 
quality control is provided in the Supplementary Information. 


Obtaining and annotating the mitochondrial genome 


To obtain the mitochondrial genome sequence, we used the iterative fine-tuning option within 
the map to reference as implemented in Geneious v.9.0.5 (www.geneious.com). This allows 
read mapping to extend past the ends of the reference sequence on each iteration. COI and 16S 
sequences of Mnais costalis (GenBank AB708347 and AF170952, respectively) were used as 
seeds against which the Illumina reads from both pools were mapped in a first step, and the 
obtained contigs were used to seed subsequent mapping, with the procedure repeated for a total 
of 425 iterations to recover the entire mitogenome. 

Mitochondrial genome annotation was done in Geneious v.9.0.5. First, open reading frames 
(ORFs) were identified, and the putative ORFs were then compared with the odonate 
mitogenomes listed above. Gene identity was confirmed by BLAST search (Altschul, Gish, 
Miller, Myers, & Lipman, 1990) against Genbank’s nr database (www.ncbi.nih.gov; date 
accessed 17 February 2016). Transfer RNA genes were identified using ARWEN (Laslett & 
Canbiick, 2008; http://mbio-serv2.mbioekol.lu.se/ARWEN/). 


Phylogenetic analyses 


Sequences of 13 protein-coding genes (PCGs), two rRNA genes and 22 tRNA genes from the 
currently available odonate mitogenomes (Table 1) were used to perform a phylogenetic analy- 
sis. Sequences from the same genes were extracted from the mitochondrial genomes available 
for Ephemera orientalis (Ephemeroptera, GenBank Accession EU 591678; Lee et al., 2009), 
Parafronurus youi (Ephemeroptera, GenBank Accession EU 349015; Zhang, Zhou, Gai, Song, 
& Zhou, 2008) and Preronarcys princeps (Plecoptera, GenBank Accession NC_006133; Stewart 
& Beckenbach, 2006); to be used as outgroups in the phylogenetic analyses. 

PCGs were aligned using Muscle (Edgar, 2004) as implemented in Geneious v.7.0.1 
(www.geneious.com). For the alignment of the 16S and 12S sequences we used Mafft v.7 
(Katoh & Standley, 2013; http://mafft.crc.jp/), with the option Q-INS-I, which considers 
the secondary structure information of the RNA. Transfer RNA sequences were aligned 
using the software LocARNA (Smith, Heyne, Richter, Will, & Backofen, 2010; Will, Joshi, 
Hofacker, Stadler, & Backofen, 2012; Will, Reiche, Hofacker, Stadler, & Backofen, 2007; 
http://mna.informatik.uni-freiburg.de/LocARNA/), which aligns tRNA sequences based on their 
sequence and structure features. Regions with ambiguities, as well as variable regions within 
each alignment, were recognised and excluded using Gblocks (Castresana, 2000; Talavera 
& Castresana, 2007; http://molevol.cmima.csic.es/castresana/Gblocks_server.html), using the 
relaxed parameters (i.e. allow smaller final blocks, allow gap position within the final blocks 
and allow less strict flanking positions). Alignments of individual genes were concatenated as 
two datasets: (1) PCG; comprising the 13 PCGs with 11,707 nucleotides; and (2) PCG+ RNA; 
comprising the 13 PCGs, two rRNAs and 22 tRNAs, with 15,110 nucleotides. Both datasets were 
used in phylogenetic analyses. 

Phylogenetic relationships were reconstructed using maximum likelihood (ML). Heuristic 
searches were carried out using the randomised accelerated maximum likelihood algorithms 
implemented in RAxML-HPC2 (Stamatakis, 2006; Stamatakis, Hoover, & Rougemont, 2008), 
through the CIPRES web portal (http://www.phylo.org). The analyses were run under the 
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GTRCAT model, and bootstrap support values were generated with a rapid bootstrap algorithm 
under the Majority Rule Criterion (autoMRE; Stamatakis et al., 2008). 


Results and discussion 


Here, we provide the third complete mitochondrial genome for a member of the family 
Calopterygidae, and further confirm the utility of next-generation sequencing to obtain whole 
odonate mitogenomes (see also Lorenzo-Carballa, Thompson, Cordero-Rivera, & Watts, 2014). 
A total of 2,813,625 short-sequence reads were mapped to obtain a single contig of 15,487 bp 
long with a mean coverage of ~ 9,000 reads. 

The size of the complete mitochondrial DNA genome of M. costalis (GenBank accession 
number KU871065) is within the size range of the other 14 complete odonate mitochondrial 
genomes (15,056—16,685 bp). It comprises the standard metazoan panel of 13 protein-coding 
genes, two rRNA genes (12S and 16S rRNA) and 22 tRNA genes (Figure 1, Table 2). Protein- 
coding genes employ the typical invertebrate mitochondrial start codons: cox1, cox2, cox3, atp6, 
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Figure 1. Genetic map of the mitochondrial genome of Mnais costalis, showing the location of the protein coding 
genes (green), transfer RNAs (pink), ribosomal RNAs (red), the intergenic spacers s1—s3 and the A +T rich region 
(grey). Picture of M. costalis by Stewart J. Plaistow. 
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Table 2. Organisation of the mitochondrial genome of the damselfly Mnais costalis. Incomplete stop codons are 
indicated with parentheses. s/—s3 are intergenic spacers. 


Feature Position Length A+ T (%) Strand Start/anti codon Stop codon 








trnI 1-72 72 65.3 + GAT - 
trnQ 68-137 70 68.5 - TTG - 
trnaM 138-206 69 66.7 + CAT - 
nad2 207-1196 990 67.9 + ATC TAA 
trnW 1194-1263 70 72.9 + TCA - 
trnC 1255-1316 62 74.2 - GCA - 
trnY 1317-1392 76 72.4 - GTA - 
sl 1393-1417 25 68 na - - 
coxl 1418-2951 1534 60.4 + ATG T(aa) 
trnL2 2951-3019 69 65.2 + TAA - 
cox2 3019-3706 688 63.2 + ATG T(aa) 
trnK 3706-3779 74 64.9 + CTT - 
trnD 3778-3842 65 78.5 + GTC = 
atp8 3843—4001 159 68.6 + ATC TAA 
atp6 3989-4672 684 64.6 + ATG TAA 
cox3 4672-5458 787 58.5 + ATG T(aa) 
trnG 5458-5522 65 76.9 + TCC -= 
nad3 5522-5875 354 67.1 + ATC TAG 
trnA 5874-5938 65 75.4 + TGC = 
trnR 5939-6005 67 61.2 + TCG = 
trnN 6003-6069 67 73.1 + GTT - 
trnS1 6070-6138 69 56.5 + GCT - 
trnE 6141-6205 65 81.5 + TTC - 
trnF 6203-6268 66 63.6 + GAA - 
s2 6269-6293 25 80 na - - 
nad5 6294-7994 1701 65.1 - ATT T(aa) 
trnH 7995-8058 64 68.7 - GTG - 
nad4 8059-9398 1339 67.9 - ATG TA(a) 
nad4L 9392-9685 294 69.2 - ATG TAA 
trnT 9687-9755 69 65.2 + TGT - 
trnP 9759-9827 69 68.1 - TGG - 
nad6 9829-10323 495 62.8 + ATC TAA 
cob 10,323-11,451 1129 6l + ATG TAG 
trnS2 11,455-11,519 65 58.5 + TGA - 
s3 11,520-11,536 17 88.2 na - - 
nadl 11,537-12,487 951 66.1 - ATG TAA 
trnLl 12,489-12,553 65 69.2 - TAG - 
l-rRNA 12,554-13,844 1291 71 - - - 
trnV 13,845-13,916 12 61.1 - TAC - 
s-rRNA 13,917-14,653 737 69.2 - = - 


A+Trich region 14,654—-15,487 834 79.7 na — — 


nad4, nad4L, cob and nadl use ATG; nad2, nad3, nad6 and atp8 use ATC, and nad5 uses ATT. 
Eight protein-coding genes have the standard stop codons TAA (nadl, nad2, nad6, nad4L, atp6 
and atp8) and TAG (nad3 and cob); coxl, cox2, cox3 and nad5 have an incomplete stop codon 
of a single T, while nad4 has an incomplete TA stop codon (Table 1). The base frequency of the 
whole mtDNA genome is A = 40%, T = 26.2%, C = 19.5% and G = 14.2%, with an overall 
A +T content of 66.2% that is within the range of other odonate mitogenomes (A +T = 64.1- 
73.1%), including the calopterygids Vestalis melania (64%) and Atrocalopteryx atrata (10%). 
All protein-coding genes have a high AT content, which varies between 69.2% (nad4L) and 
58.5% (cox3) AT content. All 22 tRNA-coding sequences of M. costalis can be folded into the 
characteristic clover-leaf secondary structure, and they range in size from 62 bp in trnC to 76 bp 
in trnY (Supplemental Information Figure S1; Table 1). 

The mtDNA genome of M. costalis is identical in gene number and gene arrangement to 
that of the other 14 complete odonate mitogenomes. Sixteen gene junctions in the mtDNA of 
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Figure 2. Maximum likelihood (ML) tree obtained when analysing nucleotides from protein coding genes, 16S and 12S and tRNAs; after ambiguous regions were excluded from all the 
alignments. Values above and below branches represent bootstrap support values for the analyses with the PCG + RNA and PCG; respectively (see main text for details). 
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M. costalis have short overlaps, with the largest junction having 13 nucleotides overlap (between 
atp8 and atp6). There are three non-coding intergenic spacers (s1—s3), which are characteristic 
of other odonate mitogenomes. The s5 spacer located between trnL and nad/ in anisopterans 
(dragonflies) is absent in M. costalis and all other zygopteran (damselflies) mitogenomes, further 
supporting the suggestion that this trait is a synapomorphy for the Zygoptera and Anisoptera (Lin 
et al., 2010; Lorenzo-Carballa et al., 2014). 

Only three SNPs (three Rs in positions 5085 [cox3]; 9638 [nad4L]) and two ambiguities (Ns 
in positions 8901 [nad4] and 14,732 [A+T rich region]) were found in the sequence, which 
indicates that some variability is expected between different individuals/morphs of M. costalis. 
Thus, having the complete mitochondrial genome for this species will help to improve primer 
design for future population genetic studies. 

Our ML analyses using both datasets (PCG and PCG + RNA) recovered the same topology. 
Contrary to the current accepted classification of the suborder Anisozygoptera as the extant sister- 
group of the Anisoptera (Dijkstra et al., 2013), our results place E. superstes as a basal group 
to the Zygoptera, although this relationship is supported with a higher bootstrap value by the 
analysis of the PCG dataset (Figure 2). These results would be in agreement with the finding that 
the mitochondrial genome of this species lacks the intergenic spacer s5 between nad/ and trnL2, 
and therefore the mitogenomic organisation of this relict odonate would be more similar to that 
of the Zygoptera (Wang ef al., 2015). 

A surprising result from our analyses is that the monophyly of the suborder Anisoptera is not 
well supported by any of the datasets (Figure 2), which could be due to the fact that the available 
mitogenomes for to anisopteran species are incomplete (Orthetrum triangulare melania and Cor- 
dulia aenea, Table 1). Obtaining more complete mitochondrial genomes for other representatives 
of the Anisoptera, as well as for more Epiophlebia species in the future, would help to improve 
the resolution of the analyses and to clarify relationships between extant odonate groups. Also, 
combining complete mitogenome data with nuclear markers would help to obtain more robust 
phylogenies for the order. 


Acknowledgements 


MOLC is funded by the EU (Marie Curie IEF PIEF-GA-2013-626504). MOLC thanks Wiebke Feindt and Kai Kamm 
for fruitful discussions and suggestions on the analyses. Two anonymous reviewers gave useful comments on an early 
version of the manuscript. 


Supplemental data 


Supplemental data for this article can be accessed at http://dx.doi.org/10.1080/13887890.2016.1234980 


References 


Altschul, S. F., Gish, W., Miller, W., Myers, E. W., & Lipman, D. J. (1990). Basic local alignment search tool. Journal 
of Molecular Biology, 215, 403-410. doi:10.1006/jmbi.1990.9999 

Ballare, E., & Ware, J. L. (2011). Dragons fly, biologists classify: an overview of molecular odonate studies, and 
our evolutionary understanding of dragonfly and damselfly (Insecta: Odonata) behavior. International Journal of 
Odonatology, 14, 137-147. doi:10.1080/13887890.2011.579538 

Cameron, S. L. (2014). Insect mitochondrial genomics: Implications for evolution and phylogeny. Annual Review of 
Entomology, 59, 95-117. doi:10.1146/annurev-ento-011613-162007 

Castresana, J. (2000). Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. 
Molecular Biology and Evolution, 17, 540-552. doi: 10.1093/oxfordjournals.molbev.a026334 

Chen, M. Y., Chaw, S. M., Wang, J. F., Villanueva, R. J. T., Nufieza, O. M., & Lin, C. P. (2015). Mitochondrial genome 
of a flashwing demoiselle Vestalis melania from the Philippine Archipelago. Mitochondrial DNA, 26, 720-721. 
doi:10.3109/19401736.2013.845757 

Córdoba-Aguilar, A. (Ed.). (2008). Dragonflies and damselflies. Model organisms for ecological and evolutionary 
research. Oxford: Oxford University Press. 

Dijkstra, K.-D. B., Bechly, G., Bybee, S. M., Dow, R. A., Dumont, H. J., Fleck, G., ... Ware, J. (2013). The classification 
and diversity of dragonflies and damselflies (Odonata). In Zhang, Z.-Q. (Ed.). Animal biodiversity: An outline of 
higher-level classification and survey of taxonomic richness. Zootaxa, 3730, 36—45. doi:10.11646/zootaxa.3703.1.9 


198 M.O. Lorenzo-Carballa et al. 


Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids 
Research, 32, 1792-1797. doi:10.1093/nar/gkh340 

Ferreira, S., Lorenzo-Carballa, M. O., Torres-Cambas, Y., Cordero-Rivera, A., Thompson, D. J., & Watts, P. 
C. (2014). New EPIC nuclear DNA sequence markers to improve the resolution of phylogeographic stud- 
ies of coenagrionids and other odonates. International Journal of Odonatology, 17, 135-147. doi:10.1080/ 
13887890.2014.950698 

Hayashi, F., Dobata, S., & Futahasi, R. (2005). Disturbed population genetics: suspected introgressive hybridization 
between two Mnais damselfly species (Odonata). Zoological Science, 22, 869-881. doi: 10.2108/zsj.22.869 

Katoh, K., & Standley, D. M. (2013). MAFFT multiple sequence alignment software Version 7: Improvements in 
Performance and Usability. Molecular Biology and Evolution, 30, 772-780. doi: 10.1093/molbev/mst010 

Kumar, C. S., Nair, R. R., Sivaramarkrishnan, K. G., Ganesh, D., Janarthanan, S., Arunachlam, M., & Sivarubam, T. 
(2012). Influence of certain forces of evolution of synonymous codon usage bias in certain species of three basal 
orders of aquatic insects. Mitochondrial DNA, 23, 447—460. doi:10.3109/19401736.2012.710203 

Laslett, D., & Canbäck B. (2008). ARWEN, a program to detect tRNA genes in metazoan mitochondrial nucleotide 
sequences. Bioinformatics, 24, 172—175. doi:10.1093/bioinformatics/btm573 

Lee, E. M., Hong, M. Y., Kim, M. I., Kim, M. J., Park, H. C., Kim, K. Y., ... Kim, I. (2009). The complete mitogenome 
sequences of the palaeopteran insects Ephemera orientalis (Ephemeroptera: Ephemeridae) and Davidius lunatus 
(Odonata: Gomphidae). Genome, 52, 810-817. doi:10.1139/g09-055 

Lin, C. P., Chen, M. Y., & Huang, J. P. (2010). The complete mitochondrial genome and phylogenomics of a damselfly, 
Euphaea formosa supports a basal Odonata within the Pterygota. Gene, 468, 20-29. doi:10.1016/j.gene.2010.08.001 

Lorenzo-Carballa, M. O., Thompson, D. J., Cordero-Rivera, A., & Watts, P. C. (2014). Next generation sequencing yields 
the complete mitochondrial genome of the scarce blue-tailed damselfly, Ischnura pumilio. Mitochondrial DNA, 25, 
247-248. doi:10.3109/19401736.2013.796518 

Plaistow, S., & Tsubaki, Y. (2000). A selective trade-off for territoriality and non-territoriality in the polymorphic dam- 
selfly Mnais costalis. Proceedings of the Royal Society of London Series B — Biological Sciences, 267, 969-975. 
doi:10.1098/rspb.2000.1098 

Simon, S., & Hadrys, H. (2013). A comparative analysis of complete mitochondrial genomes among Hexapoda. 
Molecular Phylogenetics and Evolution, 69, 393—403. doi:10.1016/J.Y mpev.2013.03.033 

Smith, C., Heyne, S., Richter, A. S., Will, S., & Backofen, R. (2010). Freiburg RNA tools: a web server integrating 
IntaRNA, ExpaRNA and LocARNA. Nucleic Acids Research, 38(Suppl.), W373-377. doi:10.1093/nar/gkq316 

Stamatakis, A. (2006). RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and 
mixed models. Bioinformatics, 22, 2688—2690. doi:10.1093/bioinformatics/btl446 

Stamatakis, A., Hoover, P., & Rougemont, J. (2008). A rapid bootstrap algorithm for the RAxML Web servers. Systematic 
Biology, 57, 758-771. doi:10.1080/10635150802429642 

Stewart, J. B., & Beckenbach, A. T. (2006). Insect mitochondrial genomics 2: the complete mitochondrial genome 
sequence of a giant stonefly, Pteronarcys princeps, asymmetric directional mutation bias, and conserved plecopteran 
A + T-region elements. Genome, 49, 815-824. doi:10.1139/g06-037 

Talavera, G., & Castresana, J. (2007). Improvement of phylogenies after removing divergent and ambiguously aligned 
blocks from protein sequence alignments. Systematic Biology, 56, 564-577. doi: 10.1080/10635 150701472164 

Tang, M., Tan, M., Meng, G., Yang, S., Su, X., Liu, S., ... Zhou, X. (2014). Multiplex sequencing of pooled mitochon- 
drial genomes — a crucial step toward biodiversity analysis using mito-metagenomics. Nucleic Acids Research, 42, 
e166. doi:10.1093/nar/gku917 

Tsubaki, Y., Hooper, R., & Siva-Jothy, M. (1997). Differences in adult and reproductive lifespan in the two male forms 
of Mnais pruinosa costalis Selys (Odonata : Calopterygidae). Researches on Population Ecology, 39, 149-155. 
doi: 10.1007/BF02765260 

Tsubaki, Y., & Okuyama, H. (2016). Adaptive loss of color polymorphism and character displacement in sympatric 
Mnais damselflies. Evolutionary Ecology, 30, 811-824. doi: 10.1007/s 10682-015-9778-3. 

Wang, J. F., Chen, M. Y., Chaw, S. M., Morii, Y., Yoshimura, M., Sota, T., & Lin, C. P. (2015). Complete mitochondrial 
genome of an enigmatic dragonfly, Epiophlebia superstes (Odonata, Epiophlebidae). Mitochondrial DNA, 26, 718- 
719. doi:10.3109/19401736.2013.845756 

Will, S., Joshi, T., Hofacker, I. L., Stadler, P. F., & Backofen, R. (2012). LocARNA-P: Accurate boundary prediction and 
improved detection of structural RNAs. RNA, 18, 900-914. doi:10.1261/rna.029041.111 

Will, S., Reiche, K., Hofacker, I. L., Stadler, P. F, & Backofen, R. (2007). Inferring non-coding RNA fami- 
lies and classes by means of genome-scale structure-based clustering. PloS Computational Biology, 3(4), e65. 
doi:10.1371/journal.pcbi.0030065 

Yamauchi, M. M., Miya, M. U., & Nishida, M. (2004). Use of a PCR-based approach for sequencing whole mitochondrial 
genomes insects: two examples (cockroach and dragonfly) based on the method developed for decapod crustaceans. 
Insect Molecular Biology, 13, 435—442. doi:10.1111/j.0962-1075.2004.00505.x 

Yu, P., Cheng, X., Ma, Y., Yu, D., & Zhang, J. (2016). The complete mitochondrial genome of Brachythemis 
contaminata (Odonata: Libellulidae). Mitochondrial DNA A, DNA Mapp. Seq. Anal., 27(3), 2272-2273. 
doi:10.3109/19401736.2014.984176 

Zhang, J., Zhou, C., Gai, Y., Song, D., & Zhou, K. (2008). The complete mitochondrial genome of Parafronu- 
rus youi (Insecta: Ephemeroptera) and phylogenetic position of the Ephemeroptera. Gene, 424, 18-24 doi:10. 
1G016/j.gene.2008.07.037. 


